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ABSTRACT 

Stellar irradiation and particles forcing strongly affect the immediate environment of extrasolar 
giant planets orbiting near their parent stars. However, it is not clear how the energy is deposited 
over the planetary atmosphere, nor how the momentum and energy spaces of the different species that 
populate the system are modified. Here, we use far-ultraviolet emission spectra from HD209458 in the 
wavelength range (1180-1710)1 to bring new insight to the composition and energetic processes in 
play in the gas nebula around the transiting planetary companion. In that frame, we consider up-to- 
date atmospheric models of the giant exoplanet where we implement non-thermal line broadening to 
simulate the impact on the transit absorption of superthermal atoms (HI, OI, and CII) populating the 
upper layers of the nebula. Our sensitivity study shows that for all existing models, a significant line 
broadening is required for OI and probably for CII lines in order to fit the observed transit absorptions. 
In that frame, we show that OI and CII are preferentially heated compared to the background gas 
with effective temperatures as large as Toi/Tb ~ 10 for OI and Tqii/Tb ~ 5 for CII. By contrast, 
the situation is much less clear for HI because several models could fit the Ly-a observations including 
either thermal HI in an atmosphere that has a dayside vertical column [HI] ~ 1.05 x 10 21 cm~ 2 , or a 
less extended thermal atmosphere but with hot HI atoms populating the upper layers of the nebula. 
If the energetic HI atoms are either of stellar origin or populations lost from the planet and energized 
in the outer layers of the nebula, our finding is that most models should converge toward one hot 
population that has an HI vertical column in the range [HI]hot ~ (2 — 4) x 10 13 cm~ 2 and an effective 
temperature in the range Tjji ~ (1 — 1.3) x 10 6 K, but with a bulk velocity that should be rather slow. 
Subject headings: line: formation — line: profiles — planetary systems — radiation mechanisms: 
non-thermal — stars: individual (HD209458) — ultraviolet: stars 



1. INTRODUCTION 

The nature of the extended nebula around giant exo- 
planets orbiting near their parent stars has remained a 
fundamental unresolved problem since it was raised a few 
years ago. Of all the exoplanets detected over the last 
decade, HD 209458b is the first discovered by the tran- 
sit technique, and the only exoplanet for which rich far- 
ultraviolet (FUV) space observations were obtained be- 
fore the Hubble Space Telescope/Space Telescope Imag- 
ing Spectrograph (HST/ST I S) ins t rument went down 
(IVidal-Madiar et al. I IpfioH l200l : IBen-.Taffel. I (pcM 
2008), hereafter, respectively, VM03, VM04, BJ07, and 
BJ08). The collected data set which covers the wave- 
length range of 1180-1710 A is a unique tool not only 
for understanding the composition and structure of the 
atmosphere of these giant planets but also for character- 
izing the processes that heat and excite them. Therefore, 
the interpretation of this data set is of great importance 
and should ultimately help answer fundamental ques- 
tions regarding the formation and evolution of extrasolar 
planets, particularly questions concerning their survival 
given their close proximity to the harsh environment of 
their parent stars. 

Historically, HD209458b was discovered transiting 
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its parent star and abs o rbing 1.5% of it s flux 
(jCharbonneau et al. I (pOOOh : iHenrv et al.. I (pOOOh ). Ex- 
tending the technique to space FUV observations at Ly- 
o (1216 A), an enormous hydrogen nebula was discov- 
ered using the HST/STIS spectrograph, covering almost 
- 8.9% ± 2.% of the stellar disc (VM03, BJ08 ). This 
corresponds to a hydrogen nebula extending over an ef- 
fective area that has at least ~ 2.47 planetary radius 
across the line of sight, tending to support the theory 
that the coronal hydrogen gas of HD209458b is inflated. 
The coincidence between the Roche lobe size across the 
planet-star line and the occultation effective area suggest 
that the gas is probably filling the exoplanetary Roche 
lobe (VM03, BJ08). Later, based on new HST/STIS 
FUV observations in the range (1180-1710A) of heavy 
elements like C and O, an occultation, respectively, of 
5. ± 2% for HI, 13 ± 4.5% for OI and 7.5 ± 3.5% for 
CII were reported (VM04). There is still a great deal 
of uncertainty regarding the strong absorptions obtained 
during transit, respectively, for the two heavy minor con- 
stituents, yet the drop-off appears comparable to the 
transit absorption obtained for HI, which, according to 
VM04, supports a blow-off scenario in the upper atmo- 
sphere of HD209458b in which an escaping element (HI) 
drags other heavy constituents outward. The problem 
is that both OI and CII are several orders of magnitude 
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less abundant than HI, and yet having the two compo- 
nents filling or even overflowing the planet's Roche lobe 
will not help to a ttain the reported transit absorptions 
([Garcia Munoz. I (|2007h ). In addition, it is not yet clear 
how wide the stellar lines are for the 01 and CII multi- 
plets because the HST/STIS observations were obtained 
in the low-resolution mod e that does not help to resolve 
the emission lines fVM04. [Garcia Munoz.1 (|2007f) V Line 
broadening processes were thus suggested to explain the 
01 and CII lines' absorption but they were never imple- 
mented in past studies. All these observational, analysis, 
and modeling uncertainties make the transit absorptions 
reported for HI, 01, and CII a real challenge for any the- 
oretical model. 

In the following, we revisit available HST archive data 
in the 1180-1710 A spectral window to conduct a sen- 
sitivity study including HI, 01, and CII abundances to- 
gether with line broadening processes in play in the up- 
per layers of the atmosphere of HD209458b. Such a 
global approach is required for any attempt to derive 
self-consistent observational constraints on the steady- 
state properties of the nebula around the exoplanet. In 
Section 2, we recall Ly-a absorption rates and profiles ob- 
tained during transit from medium spectral resolutions 
as reported by BJ08 and describe in detail the anal- 
ysis of the low-resolution observations corresponding to 
HI (1216 A ), 01 (1304A triplet), and CII (1335A mul- 
tiplet). For comparison with our model, a set of four 
constraints was selected: the absorption line profile at 
Ly-a as derived from the medium- resolution observations 
and the line-integrated absorption drop-off during tran- 
sit corresponding to HI (1216A), 01 (1300A triplet), and 
CII (1335A multiplet) as derived from the low-resolution 
observations. Our exten ded atmospheric model , based 
on the standard model of iGarcia Munoz.1 ([2007f ). is set 
forth in Section 3 and focuses on the description of the 
tidal bulge region and stellar emission lines' properties. 
In Section 4, we describe line-broadening processes that 
may affect the absorption profiles of the different com- 
ponents. Our results are discussed in Section 5 with the 
dual aims of deriving a consistent view of the gas dis- 
tribution and the processes that heat and excite the gas 
nebula around HD209458b and propounding predictions 
for future space observations. 

2. MEDIUM AND LOW RESOLUTION SPECTRAL 
OBSERVATIONS OF HD209458 

Our first data set consists of the HST/STIS archive 
observation of HD209458 obtained with the G140M 
medium-resolution grating and the 52"x0".l slit. The 
analysis of this data set was controversial, with contra- 
dictory conclusions obtained for the Ly-a transit absorp- 
tion (VM03; BJ07) before the conclusive study in BJ08 
which we adopt herein. In the global approach that we 
implement, the transit absorption spectral profile shown 
in Figure 1 is of particular interest (see also Figure 6 
in BJ08). The symmetric shape of the observed absorp- 
tion profile, due largely to the coincidence of the tran- 
sit light curves respectively from the blue and red wave- 
length regions of the Ly-a line (see Figure 3 in BJ08), led 
us to dismiss atmospheric models that produce asym- 
metric absorption profiles (see BJ08 for more details). 
In addition, considering standard atmospheric models, 



we previously showed that the hydrogen nebula around 
HD209458b should be very opaque, behaving much like 
a damped system than a classical atmosphere. For refer- 
ence, it is helpful to recall that damped systems are galac- 
tic absorbing systems that show extended Ly-a damped 
wings in the absorpti on spectra observed toward quasars 
([Wolfe et al.. I (| 1986ft ). For these systems, the [HI] col- 
umn is usually larger than 2 x 10 20 cm -2 and the tempera- 
ture is rather cool (HI thermal velocity ~ lSkms^ 1 ), sim- 
ilar to the values derived for HD209458b's nebula (BJ08). 
However, among all possible optically thick models pro- 
posed for HD209458b, no unique solution for the atmo- 
spheric hydrogen distribution could be derived as hy- 
brid models with a thin layer of superthermal hydrogen, 
may also produce satisfactory fits to the observed ab- 
sorption profile. For those reasons, in the past, we in- 
tentionally disfavored any particular model in hopes of 
obtaining additional observations in the future (BJ08). 

The second data set included in this study is the 
HST/STIS archive observation of HD209458 obtained 
with the G140L low-resolution grating and the 52"x0".2 
slit (VM04). In this mode, the HST/STIS instrument 
has a resolution clement of ~ 1.7 pixels for point sources 
around Ly-a with an average dispersion of ~ 0.6A per 
pixel. Generally, 2-3 element resolutions are assumed 
for the spectral resolution of the instrument. For an ex- 
tended source, such as the Earth's geocoronal emissions 
and the sky's background emissions which fill the whole 
slit area, the spectral resolution is ~ 8 pixels with the 
same dispersion. These details are important because 
the recorded signal in these data is a combination of a 
stellar point source emission and the geocoronal and sky 
background extended sources. The HST/STIS selected 
mode indeed shows a low-resolution spectrum but covers 
a larger wavelength band spanning a ~ 1180-1710A win- 
dow. In the present analysis, we follow exactly the same 
technique using the time-tag information as described in 
BJ08. In total, we have four visits during the planetary 
transit of HD209458, each composed of three long expo- 
sures (~ 2000 s) distributed in time over the transit pe- 
riod (see Table 1). All exposures were obtained with the 
HST/STIS G140L grating and are sampled here using a 
time bin of 300 s from the available time-tag information. 
The whole data set resulted in ~ 67 bins time series of 
the HD208458 planet-star system as a function of the or- 
bital phase angle. Coincidentally, the time sampling of 
the transit period resulted in gaps (lasting, respectively, 
~ 2493s, ~ 1169s, ~ 851s, and - 709s). The size and 
number of the gaps, along with the relatively low signal- 
to-noise ratio (S/N) of the data, make a consistent study 
of the resulting time series difficult (BJ07). Neverthe- 
less, we are able to fit synthetic light curves into each 
selected wavelength band before we derive the drop-off 
in the integrated stellar signal during transit, along with 
the appropriate statistical errors. 

For reference, a stellar spectrum is shown in Figure 
2. Aside from the emission lines previously described 
by VM04 (their Table 1), a stellar continuum must be 
subtracted before determining the transit drop-off for in- 
dividual line or band emissions (e. g., Figure 2A). For 
that purpose, we use a low-order polynomial model of 
the logarithm of the stellar continuum spectrum in or- 
der to obtain the best fit with the observed spectra in 
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TABLE 1 

HST/STIS LOW RESOLUTION DATA SET ON HD209458. ALL OBSERVATIONS WERE OBTAINED WITH TH E G140L GRATING 
AND TH E 52"x 0.2" LONG SLIT. TRANSIT CENTRAL TIME (TCT) IS DEFINED BY 2,452,826.628521 HJD (jKNUTSON ET AL.. 1 

(|2007P ). Earth's velocity is projected along line of sight in the heliocentric reference frame. 



Dataset name 


Program ID 


Date Obs. 


Time Obs. 


Start time - TCT (s) 


Duration (s) 


Earth's Velocity 


O8U601010 


10081 


2003-10-09 


03:53:41 


-11336.88 


1780. 


-15.37 


O8U601020 


10081 


2003-10-09 


05:22:40 


-5997.87 


2100. 


-15.39 


O8U601030 


10081 


2003-10-09 


05:22:40 


-237.87 


2100. 


-15.41 


O8U602010 


10081 


2003-10-19 


19:35:34 


-4436.24 


1780. 


-19.03 


O8U602020 


10081 


2003-10-19 


21:11:34 


1323.76 


2100. 


-19.05 


O8U602030 


10081 


2003-10-19 


22:47:35 


7084.78 


2100. 


-19.07 


O8U603010 


10081 


2003-11-06 


08:44:34 


-10983.48 


1780. 


-23.66 


O8U603020 


10081 


2003-11-06 


10:13:18 


-5659.48 


2100. 


-23.68 


O8U603030 


10081 


2003-11-06 


11:49:18 


100.51 


2100. 


-23.69 


O8U604010 


10081 


2003-11-24 


00:53:55 


-6709.71 


1780. 


-26.17 


O8U604020 


10081 


2003-11-24 


02:20:58 


-1486.73 


2100. 


-26.17 


O8U604030 


10081 


2003-11-24 


03:57:00 


4275.28 


2100. 


-26.17 



Velocity (km s 1 ) 
-200 -100 100 200 
1.20 b ; ; J 



X 1-10 r 
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z 0.80^ 



0.70 ^ : : 1 

121.50 121.55 121.60 
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Fig. 1. — Absorption line profile (filled circles) during transit of 
HD209458b obtained as a ratio between Ly-c* line profiles in-transit 
and unperturbed. The sky background spectral domain as defined 
here is indicated by two vertical dashed lines. The spectral window 
was restricted to ~ ±200kms — 1 from line center because the signal 
becomes much too noisy beyond. A flat absorption rate of ~ 8.4% 
fits rather well with the obtained absorption profile (da shed line). 
In BJ 08, we used the DIV1 atmospheric model of Garcia Munoz, 
( 2007) to propose a Lorentzian-like line profile model with extended 
wings that also fits quite well with these observations. 

the wavelength window ~ 1400 — 1700A (e.g., Figure 
2A). After subtraction, we find that no extra correction 
is needed for any line considered here, which lends sup- 
port to the continuum model as a good approximation 
of the stellar continuum (e.g., Figure 2B). The result- 
ing emission lines for the three selected constituents for 
the in- and out-of-transit periods are shown in Figure 3 
with the corresponding statistical errors. Stellar signal 
drop-off during transit for the emission bands are sum- 
marized in Table 2. Our set of results is consistent with 
VM04 within the large statistical uncertainties tied to 
these observations. 

In summary, the combination of HST/STIS archive ob- 
servations results in a set of constraints corresponding to 
the HI Ly-a medium-resolution absorption profile during 
transit as shown in Figure 1, and the transit absorption, 
respectively, for HI ~ 6.6% ± 2.3% in the wavelength 
band 1212-1220A, OI - 10.5% ±4.4% in the range 1299- 
1310A, and CII ~ 7.4% ± 4.7% in the range 1332-1340A 
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Fig. 2.— (A) Out-of-transit stellar spectrum of HD209458 ob- 
tained by HST/STIS G140L low-resolution grating. The time 
integration corresponds to merging exposures O8U601010 and 
O8U603010 (e. g., Table 1), similar to VM04. The stellar contin- 
uum shows the classical, solar-like near-UV signal increase trend. 
The best estimate of the continuum model of HD209458 is over- 
plotted for reference. (B) Stellar spectrum after subtraction of 
continuum. Statistical errors are shown for reference. 



(e. g., Table 2). In the following section, we will de- 
fine our atmospheric model, the different processes of line 
broadening at play in the hot atmosphere of HD209458b, 
and the shape of the stellar emission lines before we con- 
duct a simultaneous comparison to the four constraints 
assumed here. 
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TABLE 2 

Flux drop-off during transit of HD209458 derived from HST/STIS G140L low resolution observations. Typical in & out 

TRANSIT SPECTRA FOR THE DIFFERENT ELEMENTS ARE SHOWN IN FlG. 3. (a) THE SIGNAL DROP-OFF IS DERIVED FROM A BEST \ 2 FIT TO 
THE OBSERVED TIME SERIES USING SYNTHETIC CURVES FOR EACH INDIVIDUAL WAVELENGTH BAND, (b) ONE ELEMENT IS LISTED BUT 

SEVERAL EXCITED STATES MAY ALSO CONTRIBUTE. 



Atoms 






Wavelength window (nm) 


Flux drop-off a (%) 


HI 






[121.2, 122.0] 


6.6% ± 2.3% 


OI 6 






[129.9, 131.0] 


10.5% ± 4.4% 


CII 6 






[133.2, 134.0] 


7.4% ± 4.7% 


Continuum 






[140.0, 170.0] 


1.96% ± 0.42% 


-3000 


Velocity 
-2000 -1000 


(km s') 



3. 

1000 


SIMULATION OF TRANSIT ABSORPTION: 
METHOD 




o 

CO 
CD 

1 

CD 
O 

X 



0.6 
0.4 

0.2 

0.0 
-0.2 



120.5 121.0 121.5 122.0 
Wavelength (nm) 

Velocity (km s ') 
-1500 -1000 -500 500 1000 1500 




II I I I I I I I I I I I : 1IIIIHIII 



o 

CO 
CO 

1 

CD 



129.7 130.2 130.7 131.2 

Wavelength (nm) 

Velocity (km s A ) 
-1500 -1000 -500 500 1000 1500 

0.8 
0.6 
0.4 
0.2 



o 

*" 0.0 
El 

-0.2 




i I I I I I I I I I I I j -j-j I I I I I I i i I I 



132.7 133.2 133.7 

Wavelength (nm) 



134.2 



Fig. 3. — Typical in- & out-of-transit low-resolution line profiles 
respectively corresponding to HI (a), OI (b), and CII (c). Statisti- 
cal error bars are shown for the selected in-transit period. 



In this section, we take a theoretical approach in verify- 
ing whether the absorption profile for HI (e.g., Figure 1) 
and the flux drop-off during transit for the HI, OI and CII 
lines (e. g., Table 2) are predictable. To evaluate transit 
absorption spectra, integration is implemented using a 
fine grid (see Equation 1 in BJ08), taking into account 
the varying number density and temperature distribu- 
tions according to the assumed atmospheric model. We 
recall that for a target that is not spatially resolved, such 
as for transiting exoplanets, the observed signal extinc- 
tion is a weighted average of absorption from different 
regions of the exoplanetary extended atmosphere, with 
the outer layers having the maximal weight. In other 
words, because the planet is not spatially resolved, vari- 
ous sketches of gas distribution around it may lead to the 
same observed extinction (BJ08). Ingress and egress ob- 
servations, which are not yet available with high S/N, are 
required to capture more details of the gas distribution 
around the exoplanet. In our previous study (BJ08), we 
showed that asymmetric solutions in the velocity-space 
in the example of radiation pressure's induced flow out 
of the Roche lobe as proposed by VM03 or energetic neu- 
tral atoms (ENAs) produced by the interaction between 
impinging stellar w i nd p lasma and planetary neutrals 
(jHolmstrom et al~~l (|2008f l) are not compatible with the 
rather symmetric full-transit absorption as a function of 
wavelength (BJ08) if large Dopplcr-shifted absorptions 
are assumed. In the example of the ENA model, a slow, 
external, hot population could be consistent with our 
data analysis (BJ08) and will be considered here as a 
potential absorber of the Ly-a photons during transit. 
In this study, we thus focus on cylinder symmetric atmo- 
spheric models with respect to the planet-star line. Any 
departure from this ideal picture of symmetry that is not 
detected in our data set could be corrected for afterward 
starting from the solutions derived in this study. 

3.1. Atmosperic Model of Exoplanet HD209458b 

In the scenario of atmospheric inflation by stellar 
radiation and wind particles, exoplanets closely orbit- 
ing their parent stars suffer in harsh conditions that 
may lead to distortion of the out er layers of t h eir at - 
mosphere in unpredic table ways ([Yelle et al.. I (j2008); 
lLammer et al.. I (1201)9: For example, the presence of 
a stellar magnetic field may call for a magnetohydro- 
dynamics (MHD) description of its interaction with a 
possibly magnetized and out- flowing planet. Depending 
on the planetary and stellar flow regimes, the strength 
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Fig. 4. — Schematic representation of the gas distribution around 
HD209458b. The planetary gas is assumed to be bound to the 
planetary Roche lobe (regions I & II). The corresponding transit 
absorption is evaluated using the gas opacity inside that bound- 
ary. Outside that depicted region, external sources, which may be 
of stellar origin or populations lost from the planet, form region III, 
which may exist dayside, nightside, or both. Region III opacity and 
temperature are free parameters. The spherical re gion inside i 
corres ponds to the ID atmospheric model DIV1 of Garcia Munoz, 
1120071 ). In the bulge region extending from Rl"i to the limit of the 
Roche lobe (the tidal effect), density and temperature are derived 
from the DIV1 model but may be scaled independently for each 
specie. This representation sketches most potential gas distribu- 
tions in the HD209458 planet-star system, including the impact of 
the stellar wind (see Section 3.1 for more details). 

and relative orientation of the two magnetic fields, and 
the degree of ionization of the plasma, the final con- 
figuration of the interacting system could be very com- 
plex, as illustrated by similar situations in our solar sys- 
tem such as the interaction between the solar wind and 
the lo cal interstellar flow (Ratkicw icz and B en-Jaffel. 
(2002)) or the inte raction between the Jovian plasma 
and it s satellite Io (jCombi et al.. I (|1998f) : Cla rke et al.. I 
(2002)). The problem is that most models published 
thus far focus either on a one-dimensional hydrody- 
namic description that includes chemist ry but n eglect s 
the interaction with the stellar wind (lYelle. I (|2004T) : 
iGarcia MunoH (|2007t ) : iKoskinen et al. 1 (|2007l )). on hy- 
drodynamic simulations that miss most of the underlying 
micro-physics in the system but inclu de gravity and/or 
wind int eraction with the pa re nt star (Lecavelier et al.. 1 



(2004 1); iTian et al- I 112003): iSchneiter et al.' 
iMurrav-Clav et al.. I (|2009D ; iStone and Proga. 



or exclusively on the external wind pro duction of ener- 
getic neutrals (jHolmstrom et aL~l (|2008| )). 

In contrast, this study follows a semi-empirical ap- 
proach that uses a forward analysis, analyzing simple 
sketches of the gas distribution in the interaction of 
the planet-star system and determining how well they 
fit with the observational constraints. This approach is 
commonly used and has proven its robustness in several 
studies of th e upper atmosphere of g i ant planets in the 
solar system ([Gladstone et al.. I (|2004ft ; iBen-Jaffel et al.. I 
(|2007lh iSlanger et al.. I (|2008l) h Our theoretical model is 
based on sophisticated atmospheric models that include 
most of the chemistry expec ted in the HD209458 sys- 
tems (Garcia Munoz. I (|2007l )). Within that frame, we 
assume that the exoplanet is filling its Roche lobe with 
the configuration of gas distribution as shown schemati- 
cally in Figure 4. The atmosphere is elongated along the 
planet-star line by the tidal forces (this is the bulge ef- 
fect). The size of the spherically symmetric atmospheric 
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Fig. 5. — Distribution of the planetary atmosphere opacity as a 
function of the impact parameter measured across the planet-star 
line, respectively, for HI, OI, an d CII. The distribution corresponds 
to the DIV1 spherical model of Garcia Munoz, (2007) before any 
atmospheric scaling is applied. The opacity from the bulge region 
and/or from any external source should be added for the final es- 
timation of the transit absorption. 



m odel (region I) co rresponds exactly to the DIV1 model 
of lGarcia Munoz. I (|2007[ ) that assumes solar abundances. 
The main feature of this model is that it includes heavy 
elements chemistry, which directly enhances the abun- 
dance of mo st species compared to other models, as in 
lYelle. I ()2004f ) who only considered hydrogen- helium com- 
ponents. The only change we tolerate for this reference 
model is scaling the abundance of different constituents 
over the entire expanse of atmosphere. This scaling 
should encompass most models thus far published either 
for the enhancement or depletion of the abundance of the 
different species. Such abundance variation has its ori- 
gin in the chemistry thus far assumed, large-scale atmo- 
spheric processes as observed in the up per atmosphere of 
the outer pl anets of the solar system (|Ben-Jaffel et al.."1 
(fl993l [199a) . BJ08), or external sources as suggested 
by recent two-dimcntional hydrodynamic simulation of 
the planetary wind ' s inte raction with the stellar wind 
([Stone and Proga. I (|2009t l). 

To account for the tidal bulge (region II), we assume 
the DIV1 model between Kl"i (top of region I) and R^i 
(Roche lobe limit) for the number density and tempera- 
ture reference distributions. The size of the Roche lobe 
across the planet-star line should not be mu ch different 
from t h e theoretical value o f Rl"i ~ 2.89R P (jGu et al.. I 
(2003); Uaritz et alH (|2005fl ). This compares well with 
the size derived in our previous study of the stellar Ly- 
a extinction observed during transit (BJ07, BJ08). The 
DIV1 atmospheric opacities calculated for the different 
spectral bands (1216A, 1304A, and 1335A) considered 
here are shown in Figure 5. For the reference DIV1 
model, the dayside vertical column of the three species 
along the planet-star line, including the bulge region, is 
[HI] ~ 1.57 x 10 21 cm- 2 , [OI] ~ 3.94 x lO^cm- 2 , and 
[CII] ~ 2.18 x 10 17 cm- 2 . 

In addition to the thermal population described above, 
we consider two sources of hot atoms that are included in- 
dependently in our transit absorption calculations. First, 
to account for superthermal populations within the up- 
per layers of the nebula (the internal source of hot 
atoms), two free parameters are introduced, correspond- 
ing to the position of the bottom of that active region 
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in the atmosphere measured from the top of region II 
and the temperature of the hot atoms of the considered 
specie. Outside of the Roche lobe configuration sketched 
in Figure 4, we also include the possibility of an external 
source (region III) of absorbing atoms, described by its 
own atomic column and temperature. This source could 
be of stellar origin or populations lost from the planet 
dayside, nightside, or both. If the external source is of 
stellar orig in, then it coincides with t he ENA model pro- 
posed by (jHolmstrom et al.. I (|2008ft ). Along with the 
scaling factor for the whole atmosphere, this description 
should cover most scenarios of gas distribution, includ- 
ing thermal and non-thermal populations as well as those 
with or without the presence of a stellar wind. 

In the following section, in order to generate the transit 
flux drop-off, we evaluate the line profile of the different 
stellar emissions either directly from observed, unper- 
turbed stellar lines of HD209458 or from available solar 
emission lines. 

3.2. Stellar Unperturbed Emission Line Profiles 
HD209458 is GO V solar-like star (jWood et al.. I 



(|2005h ) in which a large number of emission lines orig- 
inate from the chromosphere/corona. For Ly-a , this 
is confirmed by several STIS high/medium resolution 
observations that show a self-reversal (double horn) 
shape that has ~ lA line's width, almost the same 
properties known for t he solar Ly-a line (VM03, BJ07, 
iLemaire et al.. I ((2005)). In addition, we recall that 
chromospheric Call H and K lines were observed for 
HD209458, and these are comparable to solar lines 
(jShkolnik et al- I (120051 )). If we also consider the com- 
parable size, mass, effective temperatu r e, and age of 
the tw o stars (iCharbonneau et al- I (12000ft :lHenrv et al- I 
(|200nh : IMelo et al. I (|200ffl : INordstrom et al. I (|2004ft L 
all of these factors tend to support the conclusion that 
the OI and CII solar line profiles are good estimates for 
HD209458 lines. 

In the specific case of the HI Ly-a emission line, the 
HST/STIS medium-resolution observations of HD209458 
provide a good reference for the stellar line profile (see 
Fig. 3 in BJ07). The so-called unperturbed line profile 
can be used not only for the analysis of the medium- 
resolution observations, as done in BJ07, but also can 
be convolved with the instrument line spread function of 
the HST/STIS G140L grating in the 52"x0".2 slit mode 
in order to generate a reference Ly-a line profile for the 
low-resolution data analysis. In order to verify the accu- 
racy of this method, one may start from high-resolution 
observation of HD209458 obtained with the HST/STIS 
echelle mode but at low S/N (VM03). After convolution 
with the appropriate line-spread function corresponding 
to the G140M grating, we confirm that line profiles ob- 
tained either directly or after convolution are in agree- 
ment. 

In the case of the OI and CII emission lines, the situ- 
ation is more complicated because the lines' profiles are 
not fully available. With the G140L spectra, obviously 
we have no way to access the genuine stellar lines because 
of the poor spectral resolution. Until recently, it was im- 
possible to state the origin of the OI and CII absorption 
during transit with any certainty because so much about 
the stellar line profiles that compose the OI and CII mul- 



tiplets was unknown. However, we believe the problem 
can be solved in two steps: first, we derive the weights 
of the different lines that compose each of the multi- 
plets using available HST/STIS high-resolution spectra 
of HD209458 (VM03). The data are noisy but the col- 
lected spectra allow one to make a good estimate between 
the lines' peaks in the ratio 1:1.5:1.17, respectively, for 
the OI (1302.17A, 1304.86A and 1306.03A) lines, and in 
the ratio 0.5:1.0 for the CII (1334. 53A, 1335.71A) lines. 
Note that the CII 1335. 66A line is dismissed because 
it almost merges with the 1335. 7lA line while having a 
very weak oscillator strength (e. g., Table 3). Because 
HD209458 is a solar-like star, the second step in mod- 
eling the OI and CII multiplets is to determine whether 
solar spectra of the different OI and CII lines are avail- 
able. Fortunately, OI 1302. 17A, OI 1304.53A, and OI 
1306. 7lA line profiles have bee n measured by the OSO- 
8 satellite with good accuracy (|Link et al- I JHH)). In 
addition, line profiles of, respectively, CII 1334. 53 A and 
CII 1335. 7lA were recorded by a NASA Aerobee sound- 
ing rocket (|Berger et al- I (|1970ft ) and more recently by 
the Solar and Heliospheric Observatory/Solar Ultraviolet 
Measurem ent of Emitted Radia tion (SoHO/Sumer) spec- 
trometer (|Curdt et al- I (|2001h ). The average FWHM 
of the OI and CII lines are ~ 0.2A (or ~ 45km s _1 ) 
and ~ 0.22A (or ~ 50km s" 1 ), respectively. It is im- 
portant to note that all line profiles have FWHMs far 
larger than the thermal Doppler width in the HD209458b 
atmosp here and much la r ger th an previously estimated 
(VM04; iGarcia MunolTl (|2007D ). In the following sec- 
tion, we show that our finding unquestionably confirms 
the need for an efficient line broadening in order to ex- 
plain the strong absorption observed during transit of 
HD209458. 

To obtain the stellar OI and CII line profiles at Earth's 
position, we start with the solar line profiles to which we 
apply the interstellar medium (ISM) absorption between 
HD209458 and the observer. As no observation exists 
for the local interstellar medium (LISM) OI absorption 
in t hat direction, we sca le from the column obtained for 
HI (|Wood et al.. I (12001 V using t h e \OI/ HI]tsm ~ 10~ 4 
ratio in the ISM (|Moos et al.. I (|2002fl ). For the CII 
ion column, we use the CII ion column measured to- 
ward Capela that we sc ale for the HD209458 distance 
([Wood &: Linskv.H (|1997h ). The corresponding profile of 
the unperturbed stellar Ol and CII lines at Earth's or- 
bit are shown in Figure 6. The OI and CII multiplets 
are built using the peak weights (see above) between 
the different components as derived from the HD209458 
spectrum obtained by HST/STIS in the echelle mode. 
After the lines are convolved with the appropriate line 
spread function, a comparison to the low-resolution, un- 
perturbed line profile helps, if needed, to adjust the scal- 
ing factor between the different components that make 
the multiplets. For each individual transition, a Voigt 
atomic absorption profile is assumed. Corresponding key 
parameters required to estimate atmospheric opacities 
and corresponding absorption lines are provided in Ta- 
ble 3. 
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TABLE 3 

Properties of HD209458 stellar FUV Emission Lines. These are the most important lines and continuum 

CONSIDERED IN THE PRESENT STUDY. ATOMIC PARAMETERS ARE PROVIDED AT A REFERENCE TEMPERATURE OF 10 4 K. (To IS 
THE RESONANCE CROSS SECTION AT LINE CENTER, a v IS THE VOIGT PARAMETER REQUIRED TO EVALUATE THE ATOMIC 
ABSORPTION PROFILE, AND THE PAIR (xo, Xl) ARE THE PARAMETERS THAT BETTER REPRESENT THE STELLAR LINE AS A 

double Gaussian profile ~ ^exp- +exp- J (|Link et al.. I (| 19881 )'). The parameter x is the 

photon's frequency in Doppler units, AND (xo, Xl) are provided in Doppler units at a reference temperature 

10 4 K. The double Gaussian profile only applies for those emission lines that exhibit a double horn shape 
before absorption by the ISM gas. The OI triplet is due to the transition 3s 3 5 — > 2p 3 P. The CII multiplet 

results from the transition 2s2p 2 2 D — s> 2s 2 2p 2 P°. 
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HD209458 01 Line Profiles 
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Fig. 6. — Observed solar OI (a) and CII (b) lines assumed here 
for HD209458 after absorption by the interstellar medium gas be- 
tween HD209458 and Earth. Most lines have an average width of 
~ 45kms~ 1 for the OI triplet and ~ 50kms — 1 for the CII multiplet, 
much larger than the Doppler width at the atmospheric tempera- 
ture of 1.2 X 10 4 K. 



4. SENSITIVITY OF TRANSIT ABSORPTION TO 
ATMOSPHERIC PARAMETERS AND 
PROCESSES 



4.1. Thermal Processes: Classical Approach 

In the first attempt to compare observations to the 
transit extinction predicted by our basic model, we con- 
sider thermal broadening as the only process responsi- 
ble for the formation of the atomic abs orption profiles. 
Using the DIV1 atmospheric model of iGarcia MunozTI 
(2007), the calculated transit Ly-a absorption profile is 
slightly wider than the observed profile with ax 2 ~ 1-98 
(Figure 7A), while the line- integrated transit drop-off 
at Ly-a attained by the model ~ 7.4% is also slightly 
larger than the observed drop-off (6.6% ± 2.3%), yet still 
within the measured error bars (see Table 2). The x 2 
indication is based on the medium- and low-resolution 
Ly-a observations and is provided merely to simplify the 
comparison between the multiple solutions. One way to 
decrease the atmospheric extinction is to deplete the gas 
abundance in the DIV1 model by a scaling factor over 
all of the atmosphere. For the s ake of c larity , we may 
consider the model proposed by lYelle. I (|2004D and ob- 
tain its enhancement in order to fit the observed tran- 
sit absorption. As described in BJ08, the origin of 
such enhancement in the atomic hydrogen abundance 
could b e a different chemistry that includes heavy ele- 
ments (jGarcia Munoz. I (|2007| )). HI production by col- 
lisional processes between magnetospheric particles and 
the background atmospher e similar to conditions in the 
auroral regions of Jupiter (|Ben-Jaffel et al.. 1 (|1993fl L or 
gas accumulation fro m the stellar wind that may con- 
fine the atmosphere (jStone and ProgaH fj2009h ) . In the 
example of Saturn, the origin of the HI excess is at- 
tributed to HI chemistry enhancement caused by pre- 
cipitating water by erosion from th e extended rings of 
the planet (|Ben-Jaffel et al.."1 (fl~99 5l). The atmospheric 
enhancement approach was considered in BJ08, yet it 
is only relative as one may start from an Hl-rich atmo- 
spheric model (DIV1) and consider cases of abundance 
depletion as we do here. 

Starting from the DIV1 model, the amount of deple- 
tion of the HI abundance is obtained by the applica- 
tion of scaling factor f sca that modifies the final transit 
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absorption profile at Ly-a and the corresponding signal 
drop-off. As shown in Figure 7A, the thermal broaden- 
ing of the Ly-a absorption line follows a classical curve 
of growth compared to the HI column. For reference, 
for f S ca = 1/10, 1/3, 2/3, and 1 cases shown in Figure 
7A, the corresponding \ 2 is, respectively, ~ 5.68, ~ 2.42, 
~ 1.73, and ~ 1.98. The best fit for the transit absorp- 
tion profile is obtained for a scaling factor f sca ^2/3 
with a x 2 = 1-7289. A transit drop-off of ~ 6.74% 
of the stellar flux integrated in the wavelength range 
~ 1212 — 1220A consistently fits quite well with low- 
resolution observations ~ 6.6% ± 2.3% (e. g., Table 2). 
This family of solutions based on scaling the HI abun- 
dance was initially proposed in our previous study (BJ08) 
comparing the gas nebula around HD209458b to damped 
systems. However, one of our conclusions in BJ08 em- 
phasized that no unique solution could be favored to ex- 
plain the extinction level observed in the HI Ly-a line 
profile and integrated flux. This multiplicity of solutions 
is the reason why BJ08 chose not to disseminate quan- 
titative details concerning any particular solution until 
the present in-depth study. 

The next logical step is to determine how the atmo- 
spheric model tested in the HI Ly-a analysis would work 
for the 01 triplet and the CII doublet. For that reason, 
we start from the DIV1 model using the gas configura- 
tion sketched in Figure 4 and calculate the total opac- 
ity as a function of distance from the planet's center, 
across the planet-star line. As for HI, the 01 and CII 
abundances can be independently scaled by a factor over 
the whole atmosphere. We assume that all lines that 
compose the multiplets also contribute to the absorption 
of the stellar flux. Indeed, for 01, the 1302. 17A line is 
strongly absorbed by the ISM, making it the largest con- 
tributor to the unresolved multiplet signal from the two 
remaining lines (1304.86A and 1306. 03A ) that should 
be collisionally excited (VM04). For CII, the observed 
multiplet signal should originate respectively from CII 
1334.53A and CII* 1335. 7lA the latter being collision- 
ally excited. With the opacity of 01 and CII shown in 
Figure 5, it seems quite difficult for the DIV1 model to 
fit the extinction observed for 01 and CII during transit 
even after scaling by any factor considered for HI. For 
reference, for f sca = 2/3, the line-averaged flux drop-off 
during transit for the 01 triplet is ~ 3.9%, and ~ 3.3% 
for CII. As shown in Table 2, these values fall in the 
bottom range observed for the transit absorption, mak- 
ing thermal broadening in scaled atmospheric models a 
marginal solution for the 01 and probably CII lines. 

The challenge we face with 01 and CII transit absorp- 
tions is twofold. Indeed, assuming the transit absorption 
shown in Table 2, and despite the large statistical un- 
certainties, how is it that minor constituents that are 4 
orders of magnitude less abundant than HI show similar, 
if not stronger, absorption during transit? In addition, 
the transit absorption seems even stronger for neutral 01 
than for ion CII. It follows that for the 01 and CII ab- 
sorption lines in the HD209458b extended atmosphere, 
another explanation is required that may compensate 
for their weaker abundances compared to HI and that 
might then provide comparable absorption drop-off dur- 
ing transit for the three constituents. The missing mech- 
anism that may explain the high absorption observed for 
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Fig. 7. — (a) Comparison of observed Ly-c? absorption line pro- 
file (filled circles) during transit of HD209458b to theoretical mod- 
els assuming only thermal broadening as a function of the atmo- 
spheric scaling factor fsca- (Dotted line) f sca = 1, (solid line) 
fsca = 2/3 represents the best fit in terms of x 2 f° r these kinds 
of models, (dashed line) f sca = 1/3, and (dotted-dashed line) 
fsca = 1/10. (b) Comparison between observed (filled circles) 
and theoretical best fit Ly-ct absorption line profile during tran- 
sit of HD209458b obtained for: (solid line) atmosphere scaled by 
fsca = 2/3 with thermal broadening alone, (dotted line) atmo- 
sphere scaled by f SC a = 1/3 with energetic HI (Thi/Tb = 3.5) 
in region II, and (dashed line) atmosphere scaled by f SC a = 1/10 
with energetic HI confined to region III ([HI]hot = 4. X 10 13 cm — 2 
and Thi/Tb = 102). Note the spectral inflection that appears 
at ~ ±100km/s in the wings of the Ly-a absorption profile when 
superthermal HI atoms are included. 

both OI and CII lines may be related to broadening pro- 
cesses that should be in play in the upper atmosphere of 
HD209458b. Thi s was d iscussed qualitatively in the past 
(|Garcia Munoz. I ([2007)) but no attempt was made to 
evaluate the spectral broadening that may explain the 
large absorption observed during transit. In addition, 
the unknown spectral shape of the stellar lines at se- 
lected wavelengths was a source of great uncertainties on 
modeling the transit absorption of OI and CII. In the 
following, we elicit key information on the shape of the 
stellar lines and the unconstrained broadening processes. 

4.2. Non-thermal Processes: Global approach 

Here we consider two distinct sources of energetic 
atoms, depending on the region where the particles are 
present. For instance, we assume that all energetic atoms 
produced inside the Roche lobe limit form the so-called 
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internal source, while those produced outside represent 
the external source. In the frame of our description of 
the exoplanet's atmosphere (e.g., Figure 4), the internal 
source could be present in regions I and II, while the 
external source is confined to region III. 

Generally, the existence of energetic particles results 
from the departure of their velocity distribution function 
from equilibrium due to the imbalance between pe rturb- 
ing and restoring processes (jMever-Vernet. I (|200ll )). For 
neutrals, the sources of energetic atoms are chemical- 
and radiation-induced reactions and/or energization by 
atmospheric meso-turbulence. For partially ionized 
plasma, the mechanisms of energization could origi- 
nate from wave-particle interactions or mesoscale turbu- 
lence. The resulting velocity distribution varies from one 
constituent to another, greatly enhanc ing the br o aden- 
ing e ffect for particles of larger mass (jSchizgal. I (HH, 
120071) ). For reference, this mass dependence corresponds 
to the so-called preferential ion heating and acceleration 
with respect to hydrog en in the extended solar corona 
(|Cranmer et al.. I ((2008)) . The problem is complex and 
requires solving Fokker-Planck equations for the parti- 
cles' velocity distribution function that is sensi tive to the 
neutral or ch arged state of th e background gas (Schizgal, 
(|2004 120071) . For example, iSchizgal. I ((20071 ) evaluated 
the mass dependence of the velocity distribution function 
to wave-particles' interactions that may be quite different 
from a Maxwellian or even a kappa distribution. Meso- 
turbulence may also broaden absorption lines in a very 
complex way, particularly when the gas departs from a 
steady state (jLoucif fc Ben-Jaffel. I (|2002l) ). 

Our intent is not to conduct such heavy modeling, but 
rather to follow the rich heritage of the many studies of 
the at mosphere of solar system's planets and th e helio - 
sphere (ISchizgal fc ArkosTl (11996l):IShematovich~l (2004); 
iFahr et al- I (|2007[ k iNagv. I (|2008t )T TZank et al. I $00$) . 
and references therein). Energetic atoms are commonly 
detected in these objects at all distances from the Sun. 
When we recall the huge stellar energy available at the 
orbit of giant explanets, we are not persuaded that a 
sound analysis can neglect or ignore such energetic atoms 
in their extended atmospheres. Similar to prior s tudies 
conducted on the solar coron a (TCranmer et al., I (119981 
200 8])) or the J o vian corona HBen-Jaffel et al- I (11993 



Emerich et alTI (|1996t) ; iGladstone et al- I (|2004T )V the 



preferential broadening for the OI and CII lines is ac- 
counted for by assuming synthetic velocity distributions 
with one parameter corresponding to the effective tem- 
perature T{ of the species i that should be derived from 
the comparison to the transit absorption observed for the 
selected constituent. 

4.2.1. Impact of internal Source of Energetic Atoms 

In our globa l approach, we s t art fr om the DIV1 stan- 
dard model of IGarcia Munoz.l (|2007t ). but adopt values 
of Uca (1, 2/3, 1/2, 1/3, 1/10) for scaling OI and CII 
abundances in the atmosphere. We do not consider larger 
values of f sca as models with f sca ~ 2/3 — 1.0 already 
correspond to an enriched atmosphere and fit rather well 
with the transit absorption at HI Ly-a (see Figure 7A). 
Independent scaling of each species may help to locate 
any departure from solar mixing ratios that are assumed 
deep in the atmosphere. 

The region where collisional and wave-particles' broad- 



ening is acting is placed on top of the atmosphere and 
defined by the bottom position, the first free parame- 
ter. In such an active region, the velocity distribution of 
an individual constituent is allowed to depart from the 
background. The departure from equilibrium is defined 
by the ratio Ti/Ts of the specie's effective temperature 
to the background temperature. The results of the com- 
parison of our model to the observed transit absorption 
are provided as a function of the atmospheric scaling fac- 
tor adopted and the position of the bottom of the active 
layer. Our results are summarized in Table 4 and are dis- 
cussed below to derive key properties of hot populations 
produced in regions I and II. 

First, a quick survey of Table 4 shows that a strong dif- 
ferential heating of OI and CII is required with respect to 
the background atmosphere in order to fit the assumed 
OI and CII transit absorptions. This result is certain be- 
cause preferential atomic heating is required for OI and 
CII in all possible atmospheric cases in order to recover 
the transit absorptions observed. In addition, the derived 
differential heating is here found to be stronger for OI 
than for CII, which nicely follows the mass depend ence 
as derived by theoretical stud ies | Schizgal.l (120071) an d 
observed in the solar corona ((Cranmer et al- I (|2008l) ). 
However, the corresponding absorption line broadening 
depends on the assumed thickness of the non-thermal 
region and on the scaling factor of the abundance of 
the species considered over the whole atmosphere. Ob- 
viously, our set of observations cannot unambiguously 
discriminate between the different solutions, yet all of 
them show a need for a preferential heating of OI and 
CII compared to the background atmosphere (e. g., Ta- 
ble 4). 

To illustrate the degeneracy of solutions, we consider 
the HI (1215.67A), OI (1304.86A), and CII (1334.53A) 
absorption profiles used to fit the transit absorption with 
our standard atmosphere scaled by a factor f sca ~ 2/3. 
For HI, no internal hot hydrogen is required because ther- 
mal broadening is enough to bring the model to the de- 
sired transit absorption at Ly-a (e.g., Figure 7A). Assum- 
ing a hot layer confined to a region above r ~ 213 x 10 3 
km (~ 2.25 Rp), a differential heating T i/T B ~ 10.75 is 
required for Ol in order to fit the transit absorption ob- 
served around ~ 1304A, while for CII, T C ii/T b ~ 5.5 is 
needed to recover the transit absorption assumed around 
~ 1334A. Now, for depleted atmospheres corresponding 
to smaller scaling factors and for the same position of the 
active layer, more efficient preferential broadening is re- 
quired for OI and CII lines to balance the loss in absorp- 
tion due to smaller atomic opacity. As shown in Figure 8, 
for an atmospheric scaling factor ranging from 2/3 down 
to 1/10, Toi/Tb varies from ~ 10.75 to ~ 32 for oxygen, 
while for CII, Tcii/Tb increases smoothly from ~ 5.5 up 
to ~ 12.25. In other words, the smaller the atmospheric 
scaling, the stronger the non-thermal broadening needed 
to fit the observed transit absorptions. Intuitively, this 
result simply reflects the need for an extra line opac- 
ity to compensate for reduced atomic opacity, an extra 
opacity that is produced by a larger effective tempera- 
ture. For HI, the fit obtained with an internal source 
of hot HI is rather poor compared to the fit obtained 
for thermal broadening (see Figure 7B). Moreover, the 
corresponding Ly-a spectral shape has a pronounced U- 
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TABLE 4 

Summary of sensitivity study versus atmospheric scaling factor and position of hot atoms sources required to fit 
HST/STIS archival FUV observations (see section 4). Cases for which thermal and non-thermal line broadening could not 

EXPLAIN ASSUMED OBSERVATIONS ARE INDICATED BY "NO" FOR THE CORRESPONDING SPECIE. HOT ATOMS SOURCES ARE CONFINED ABOVE 
THE INDICATED R 3 POSITION, WHERE R„ IS THE PLANETARY RADIUS. THE FIRST FIVE POSITIONS ARE FOR INTERNAL SOURCES (REGION I AND 

II). WHILE THE LAST IS FOR. AN EXTERNAL SOURCE (REGION III). 
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shape that seems different from observations with x 2 > 2 
in all cases, making this process an improbable explana- 
tion for the observed transit absorption at Ly-a . Note, 
however, that if the internal hot atoms are placed high 
in the atmosphere, close to the Roche lobe edge, then it 
would be difficult to distinguish them from the external 
source, a scenario discussed below. In summary, while 
a standard model with ^2/3 abundances gives the best 
X 2 fit for HI, along with a reasonable picture of the dif- 
ferential line broadening between HI, OI, and CII (e.g., 
Figure 8), other models with depleted OI and CII abun- 
dances (compared to the reference DIV1 model) provide 
satisfactory fits to observations but with rather stronger 
broadening needed for OI and CII absorption lines (see 
Figs. 8 and 9). Later, we show how an external source 
of hot atoms may modify this picture for HI but not for 
OI and CII. 

After the atmospheric scaling factor, the second impor- 
tant parameter is the size of the atmospheric active layer, 
where non-thermal line broadening is operating preferen- 
tially on OI and CII lines. To explore the sensitivity of 
the effective temperature to the thickness of the active 
layer (or, equivalently, to the position R s of its bottom), 
we vary the latter position above and below Rl"i of the 
bottom boundary of the bulge region (e.g., Figure 4 for 
the atmospheric structure). Assuming our standard at- 



mospheric model with f sca ~ 2/3, the transit absorp- 
tion profiles obtained for an active layer positioned above 
R s = 213 x 10 3 km (~ 2.25R P ), 243 x 10 3 km (~ 2.57i? p ) 
, and 273 x 10 3 km (~ 2.89i? p ) are shown in Figures 9A 
and 9B, respectively, for OI and CII lines. For reference, 
a line profile with the only case of thermal broadening 
is also shown. It is interesting to note distinct features 
in the wings of the line profiles due to the change of the 
position R s . Indeed, when hot atoms are confined high 
in the atmosphere, the absorption profile bears two dis- 
tinct signatures, one narrow distribution from the ther- 
mal component around the line core spectral range and 
a second, more extended distribution in the far wings 
that originates from hot atoms. It is interesting to note 
that the jump from one distribution to another appears 
as a inflection in the absorption line profile that may be 
confused with a Doppler shifted spectral signature from a 
putative fast moving population around ~ ±100km/s for 
HI and ~ ±20km/s for OI and CII. The extra-spectral 
broadening induced by the change in the hot source po- 
sition also translates into a shift of the effective tem- 
perature as shown in Figure 8 (dashed lines). In any 
case, the position of the active layer for either OI or CII 
could not be very high; otherwise, the required heating 
would be excessive (e. g., Table 4). From another side, 
it is difficult to maintain a hot layer very deep in the 
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Fig. 8. — Ratio of effective temperature to background gas tem- 
perature (Tb) for the two species considered, OI and CII, as a 
function of the atmospheric abundance scaling factor. The scal- 
ing factor is applied to the whole atmosphere taken as the DIV1 
model from Garcia Munoz, (2007). (Solid line) Energetic atoms 
confined to region above 213 X 10 3 km (~ 2.25 R P ), a position al- 
most 6 X 10 4 km below R^i/j. (Dashed line) same but for a region 
confined to a position 3 X 10 4 km above Ryi i . Changing the bottom 
of the hot atoms region not only shifts the whole effective tempera- 
ture distribution but also may make the region so optically thin for 
particular values of f SC a that no effective temperature could pro- 
duce the line broadening required to fit the transit absorption. In 
the example shown, OI has no hot source solution for f sca = 1/10 
and a source position above ~ 2.57R p . 



atmosphere because restoring processes become impor- 
tant with increasing pressures, making it unlikely that 
any forcing process can keep a large differential heating 
between constituents in a stratified atmosphere. 

4.2.2. Impact of External Source of Energetic atoms 

The last parameter of importance in our study is the 
impact of an external hot atomic source. First, we recall 
that the spectral signature of such a population should 
show little to no Doppler shift to be consistent with the 
medium-resolution Ly-a observation (BJ08). Second, we 
recall our initial conclusion that the atmospheric scal- 
ing factor should not exceed 2/3, a value for which no 
hot population is required to recover the observed tran- 
sit absorptions at Ly-a (e. g., Section 4.1). Third, the 
origin of the external source could be stellar or popula- 
tions lost from the planet over time. With that in mind, 
five values of the atmospheric scaling factor f sca are con- 
sidered in the range 1/10-2/3. Then, for each value of 
fsca, we vary the thickness and temperature of the exter- 
nal hot layer as free parameters in order to minimize the 
X 2 corresponding to the Ly-a observations. For all val- 
ues of fsca, our study shows that having an additional, 
external layer of hot HI atoms yields transit absorptions 
that compare well with flux-integrated and line profile 
Ly-a observations but with rather different spectral sig- 
natures that could not be constrained from the only data 
available. As shown in Figure 7B, the external source, 
the internal source, and the source-free cases show dis- 
tinct spectral signatures in the absorption profile during 
transit that should help to discriminate between the dif- 
ferent solutions when high-resolution and better S/N ob- 
servations become available. Despite this degeneracy, our 
study shows that if the atmospheric scaling factor f sca is 
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Fig. 9. — (a) Model of OI (130.486 nm) absorption profile during 
full transit estimated for different positions of the hot OI atomic 
source for an atmosphere scaled from the DIV1 model by a factor 
of 1/3. For reference, the solid line shows the transit absorption 
profile obtained with the same atmospheric model but without 
the inclusion of hot atomic OI. The dashed line shows a hot OI 
source confined above ~ 2.25 R p and an effective temperature of 
Toi/Tb ~ 13.75; the dotted line shows a hot OI source confined 
above ~ 2.57R P and an effective temperature of Toi/Tb ~ 20.75; 
and the dotted-dashed line shows a hot OI source confined above 
~ 2.83-Rp and an effective temperature of Tqi/Tb ~ 76. (b) Model 
of CII (133.453 nm) absorption profile during full transit estimated 
for different positions of the hot CII atomic source for an atmo- 
sphere scaled from the DIV1 model by a factor of 2/3. For ref- 
erence, the solid line shows the transit absorption profile obtained 
with the same atmospheric model but with or without the inclusion 
of hot atomic CII. The dotted line shows a hot CII source confined 
above ~ 2.57R P and an effective temperature of Tqii/Tb ~ 7, 
(dashed line) Hot CII source confined above ~ 2.89-Rp and an ef- 
fective temperature of Tcn/Tg ~ 16, and (dotted-dashed line) 
Hot CII source confined above ~ 3.21iJ p and an effective tempera- 
ture of T C ii /T B ~ 28. Note that for both OI and CII absorption 
profiles, a spectral inflection appears at ~ ±20km/s in the wings 
when hot atoms are included. 



strictly smaller than 2/3, a population of hot HI atoms is 
then required. If this population is confined on top of the 
Roche lobe limit, its parameters (hot HI column and ef- 
fective temperature) should follow the distribution shown 
in Figure 10. Interestingly, for most solutions, the effec- 
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Fig. 10. — External hot HI source (region III) parameters' 
([-ff/]hot! Thi/Tb) sensitivity to atmospheric scaling factor of 
(fsca)- For each value of fsca, the hot HI layer parameters are 
derived from the best fit to Ly-a observations in terms of \ 2 ■ F° r 
fsca = 2/3, no hot HI atoms are required. 

tive temperature seems to fall around Thi /Tb ~ 90—110 
(or T HI ~ 1.2 x 10 6 K for T B = 1.2 x 10 4 K) while the hot 
HI column is in the range [HI]hot ~ (2 — 4) x 10 13 cm~ 2 . 
It is remarkable to see that, if assumed of stellar origin, 
the external HI source derived here compares nicely wit h 
the ENA solution proposed bv lHolmstrom et al.. I (|2008[) . 
the only difference being that their bulk speed should be 
weak. 

By contrast, for OI and CII lines, our results show 
that the transit absorption is insensitive to the presence 
of hot atoms in region III. This non-sensitivity could be 
explained by the fact that region III is optically trans- 
parent for all UV lines of heavy elements considered here 
because the corresponding abundances are greatly de- 
pleted near the Roche lobe limit (~ 10 3 — 10 4 cm -3 ). For 
these reasons, such an external population, if it exists, 
could not on its own explain the strong absorption ob- 
served for OI and CII unless an additional non-thermal 
line broadening is included yet is deeper inside regions I 
and II for the two constituents (see the previous section). 

5. EXOPLANET-STAR ENVIRONMENT: ROLE OF 
ENERGETIC ATOMS 

At this level, it may be worth casting our result in the 
general context of a giant exoplanet's interaction with the 
harsh environment caused by its nearby parent star. In 
that context, the roles of the stellar wind and its related 
magnetic field are not clear, particularly in regard to 
their merging with the planetary flow and magnetic field. 
Assuming that hot planets are locked with the same side 
facing the star, some studies derived that, contrary to 
our quickly-rotating Jupiter, the strength of their mag- 
neti c field should be very weak, or it may no t exist at 
all (lip et al. I (12001: ISanchez-Layega. I (|200l : iZarka. I 
(2007): iGriegmeier. J.-M. et al.. I (j200% ). In that con- 
text, hot atomic populations should be produced in the 
system with an energy spectrum and a spatial distribu- 
tion that should both depend critically on the assumed 
magnetic configuration in the plasma. As illustrated by 
comparable situations in our solar system such as the in- 
teraction between the solar wind and the LISM, or the 
interaction between the Jovian magnetospheric plasma 



and the outflow from the satellite Io, one may legiti- 
mately postulate that detection of hot atoms from an 
exoplanetary system should reveal new pieces of infor- 
mation regarding the spatial distribution of the plas- 
mas and magnetic fields prese nt and the strengt h of th e 
interaction processes in play (IBen-Jaffel et al.. I |2000); 
Ratkiewicz and Ben-JaffeTI (12002ft : iZank et al.. I (2009); 



Combi et al.. 1 1 1199a ). While the diagnostic will not be 



easy, we believe it is fundamental for deriving key pa- 
rameters of the star-planet magnetic environment. It 
follows, then, that a self-consistent modeling of the inter- 
acting planet-star system, though difficult, will be much- 
needed in the near future. The final configuration will be 
complex and will surely require three-dimensional MHD 
modeling of the interacting plasmas. 

In our most recent study, we discussed most of the ex- 
isting models proposed to describe the gas distribution in 
the HD209458b system based on HST Ly-a observations 
(BJ08). More recently, one-dimensional and two- 
dimensional hydrodynamic simulations have been pro- 
posed in the literature and draw the interesting conclu- 
sion that the stellar wind may, under specific conditions, 
confine the planetary flow and ther eby possibly enhance 
the to t al opacity during transi t (|Murrav-Clav et aT~l 
pOOl ; iStone and Proga. I pOOl ). This possibility of- 
fers an additional, intriguing process for enhancing the 
system opacity from outside compared to the internal 
enhancemen t induced by he a vy ele ments chemistry as 
described in IGarcia Munoz.l (|2007D . Such a possibility 
could prove to be of no small importance in uncovering 
the balance between thermal and non-thermal processes 
in shaping the transit absorption for all species, particu- 
larly for OI and CII. 

The problem with most models proposed in the lit- 
erature is that they fail to recover the observed Ly- 
a transit absorption. For instance, starting from a one- 
dimensional hydrodynamic description of the interaction 
of a planetary wind wit h a stellar wind but missing the 
atmospheric chemistry, iMurrav-Clav et al.. I (|2009l ) de- 
rived a transit absorption, corresponding to the wave- 
length range defined by VM03 and BJ08, of - 2.9% 
that falls short of the 8.4% — 8.9% absorption reported 
in BJ08. These authors thus questioned the accuracy of 
the transit absorption derived at this point from the HST 
observations, arguing that the transit Ly-a absorption re- 
ported from other low-resolution observations was much 
weaker (VM04). Modeling the fully- integrated Ly-a line 
absorption during transit, they thus predicted ~ 2.4% 
absorption that still falls short of the absorption level 
reported thus far. With in that frame of reasoning, 
IMurrav-Clav et al.. I (|2009f ) thus concluded that the ori- 
gin of the missing absorption in their wind model at 
ilOOkms -1 from line center was uncertain and that a 
good candidate could be the charge exchange between 
the planetary and stellar winds. However, we recall that 
an absorption at the "magic number" of ±100kms _1 , re- 
peatedly discussed in several studies since the first report 
of the HST observations in 2003, could be the signature 
of a population flowing at that velocity range, yet it could 
also be an opacity effect purely related to a simple curve 
of growth broadening of an absorption line by the gas 
opacity (BJ08). For reference, the wavelength range re- 
ported in VM03 and BJ08 was virtually dictated by the 
wavelength range where the HST spectrum has a good 



Ben-Jaffel & Hosseini 



13 



S/N and is not truly directly related to the observed ob- 
ject. Secondly, as we show here (see Section 3), the low- 
and medium-resolution observations are now compatible 
and there is no opposition between the transit absorption 
line profile shown in Figure 1 and the transit absorption 
reported for the integrated line (e. g . j Tab le 2). It fol- 
lows that if the iMurrav-Clav et aTTTT ([2009D predictions 
for the transit Ly-a absorption fall short of the obser- 
vations, respectively, from medium- and low-resolutions, 
this only means that their model does not render an ac- 
curate picture of the gas distribution in the HD209458 
system. As a matter of fact, when scaling the hydrogen 
content in their model by a factor of 3 — 5 as proposed 
by BJ08, these authors could obtain a transit absorption 
comparable to observations. However, invoking charge 
exchange without any restriction is probably not a good 
idea if the stellar wind speed is strong enough to leave a 
spectral signature that is not observed (BJ08). 

In contrast to models previously described, here we 
take a forward analysis approach that proves very use- 
ful, at least for translating the transit absorptions ob- 
served for different species in the HD209458 system into 
atmospheric parameters and processes. First, our anal- 
ysis of the FUV HST/STIS archives confirms relatively 
strong absorptions previously reported for 01 and CII 
but with larger statistical errors. For the heavy ele- 
ments, our analysis shows that a preferential heating 
is required, which is consistent with the energy balance 
in play due to the proximity of the parent star, either 
through tidal disturbances or the impact of the stellar 
radiation field and the stellar wind particles and mag- 
netic field. Indeed, hot atoms are minor populations that 
receive energy from both the atmospheric pool of the 
HI dominant species and the stellar impinging particles 
and photons. The two sources are connected but have 
enough energy to over-heat minor components thanks 
to the imbalance between per turbing and restoring pro- 
cesses (|Mever-Vernet. I (|200lD V For example, assuming 
only 1% of the thermal energy stored in the background 
gas (at the background temperature Tb ~ 10 4 K) is 
transmitted via non-equilibrium processes to minor con- 
stituents (mixing ratio ~ 10~ 4 ), it is not difficult to find 
that minor elements could be over-heated by large factors 
that may explain the temperatures reported here. In ad- 
dition, the large line broadening required to recover the 
observed transit absorptions in all lines of heavy species is 
consistent with the theoretical prediction of preferential 
heating of heavy species in a non-equilibrium gas. For 
01, the ~ 2.3(7 transit absorption detection requires hot 
01 with an effective temperature ratio of Toi/Tb ~ 10 
or more, depending on the thickness of the active region. 
For CII, the transit absorption comes with larger sta- 
tistical errors, yet it also requires hot CII ions with an 
effective temperature ratio of Ten /Tb ~ 5 or more, but 
definitely smaller than the temperature ratio obtained 
for 01. As shown in Figure 8, the ratio between Toi and 
Ten seems consistent with theoretical predictions about 
the mass dependence of the broadening velocit y distribu- 
tion o f heavy elements in disturbed plasmas (jSchizgal. I 
(|2007| ) ) . Conversely, starting from such theoretical mod- 
els, if non-thermal heating is detected for the 01 atoms 
(with a mass M = 16), as here, then we must expect 
similar heating of the CII ions (with a mass M = 12), 
but not as strong as that obtained for 01. With that 



in mind, we may confidently state that both observa- 
tional and theoretical results tend to support a scenario 
of preferential heating of 01 and CII with respect to 
the background atmosphere of HD209458b. The answer 
to exactly how much and where atomic energization (or 
non-thermal heating) is operating will require future ob- 
servations with better S/N and spectral resolution. 

For the dominant component HI, the observations are 
better quality, yet the conclusions are still uncertain. 
Our sensitivity study favors two families of solutions for 
HI (see section 4). The first family requires only thermal 
HI in a thick atmosphere with a large dayside vertical 
column, while the second family is based on Hi-depleted 
atmospheres with respect to the DIV1 model (depletion 
larger than 3/2), but includes a thin source of hot HI 
confined to the outer layers of the Roche lobe limit of 
the nebula. If the second solution is confirmed (see Fig- 
ure 7), high-resolution and S/N Ly-a absorption profiles 
observed from ingress to egress orbital positions should 
provide the related orbital distribution of energetic neu- 
tral HI (effective temperature and HI column), which 
in turn should help map the spatial distribution of the 
energization processes in play in the system. In addi- 
tion, if the hot HI source is of stellar origin, our results 
place strong constraints on the stellar wind and should 
require it to be rather slow, consistent with the picture 
of the exoplanet evolving in the stellar extended corona 
where t he flow still may be slow compared to farther dis- 
tances (jPreusse et al.. I (|2005t ): lGrie/3meier. J.-M. et al. I 
(2007)). Finally, we note that even if the external source 
appears to be a potential option for HI, it is not a viable 
solution for the less-abundant heavy species 01 and CII, 
for which an energization process is required deeper in 
the atmosphere. 

At this stage, one may wonder what diagnostic can 
assist in obtaining the thickness and effective tempera- 
ture of the regions where the different energetic atoms 
are confined. First, we note that the preferential heating 
described here should not be uniformly distributed in a 
confined region, but rather should be distributed with al- 
titude. Now, depending on the depth at which hot atoms 
are present, our model calculations seem to indicate that 
all HI, 01, and CII absorption line profiles predicted dur- 
ing full transit should show distinct shapes, particularly 
in the line wings (see Figures 7 and 9). This means that 
high-resolution and S/N observations obtained for dis- 
tinct species (the mass effect) may help discriminate be- 
tween the different solutions. Of particular interest, the 
01 and CII transit absorption profiles should not only 
constrain the balance between thermal and non-thermal 
atomic populations for these species, but also help to 
answer the question of whether heavy elements in hot 
Jupiters are depleted with respect to solar abundance 
values. In addition, ingress and egress observations with 
high S/N should provide more details about the spatial 
distribution of the occulting plasma, which in turn should 
bring new insight as to the magnetic configuration and 
energization processes via the inferred abundance and 
distribution of the energetic populations. 

6. SUMMARY AND CONCLUSIONS 

HST archive low- and medium-resolution observations 
of the HD209458 system in the FUV spectral domain 
are revisited to uncover the gas distribution and ki- 
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netic processes in the nebula that envelops the exoplanet. 
Within that framework, three important constituents- 
HI, 01, and CH-are considered. Our analysis of the 
low-resolution observations obtained in the wavelength 
range of 1180 — 1710 A shows a transit absorption of 
~ 6.6% ± 2.3% for HI, - 10.5% ± 4.4% for 01, and 
- 7.4% ± 4.7% for CII. Past results from the anal- 
ysis of medium-resolution observations for the HI Ly- 
a component are utilized here (BJ08). Starting from 
a standard model that includes most of the important 
chemistry in play in an atmosphere forced by radiation, 
particles, and gravity from the star, we built a more so- 
phisticated model to account for broadening processes 
that may originate from superthermal atoms in the upper 
layers of the nebula or that may be of stellar origin. To 
estimate the transit effect, the stellar line profiles must 
be known. While the stellar Ly-a line profile is known to 
a high degree of accuracy from HST medium- and high- 
spectral resolution observations, such is not the case for 
the 01 and CII lines for which existing solar line profiles 
are assumed. The major difference between this study 
and previous ones is consideration of 01 and CII lines 
simultaneously with the HI Ly-a line. The challenge was 
to obtain for 01 and CII almost the same level of tran- 
sit absorption as for HI when the abundance is at least 
3 or more orders of magnitude smaller for the heavy el- 
ements, particularly when the atomic absorption cross- 
sections for the three sets of lines have the same order 
of magnitude. A differential process that must compen- 
sate for the OI/HI and CII/HI lower abundances is then 
required in order to determine the comparable magni- 
tude of transit absorption by the three species. Here, we 
propose particles-particles collisions, wave-particles in- 
teraction, and mesoturbulence effects that may enhance 
line broadening during the stellar radiation transmission. 
Such processes are driven by huge forces caused by the 
impact of stellar radiation, impinging particles, and tidal 
and magnetic distortions, all of which act cumulatively 
on the exoplanetary atmosphere. 

Three families of models are thus considered: first, a 
case free of energetic atoms; second, a case with hot 
atoms included internally (inside regions 1 and 2); and 
third, a case with external hot atoms (within region III) . 
With an independent scaling factor for each specie, this 
description covers most models discussed or proposed to 
date in the literature. First, the comparison of our the- 
oretical model to HST observations of HD209458 shows 
that for all models considered, a population of hot 01 and 
CII atoms, confined inside the Roche lobe, is required in 
order to fit with the transit absorption observed for the 
two species' lines. Our parametric analysis clearly shows 
that for all considered cases, a preferential heating of 01 
and CII compared to HI is required, with a magnitude 
that depends on the assumed atmospheric model (e.g., 
Figures 8 and 9) but with an effective temperature that is 
definitely higher than the background temperature (e.g., 
Figure 8) . In addition, all effective temperatures derived 
here are higher for 01 than for CII, consistent with the 
theoretical prediction of the mass dependence of veloc- 
ity distr ibution of heavy elements in a non-equilibrium 
plasma (jSchizgal. I ■ 20071)). Interestingly, the presence 
of energetic atoms may translate into a spectral inflec- 
tion that appears ~ ±20km/s from the line center of 



the 01 and CII lines (Figure 9). Future detection of 
this spectral feature would confirm the direct signature 
of superthermal atoms populating the upper atmosphere 
ofHD209458b. 

For HI, by contrast, two families of models provide a 
satisfactory fit to available Ly-a observations (see Fig- 
ure 7). For instance, when considering a model free of 
energetic HI atoms, it is possible to fit both the tran- 
sit absorption profile and the integrated flux drop-off at 
Ly-a simultaneously with an atmosphere that has a day- 
side column [HI] ~ 1.05 x 10 21 cm~ 2 , corresponding to 
a depletion in HI by a factor 3/2 wi th respect to the 
DIV1 model of iGarcia Munoz. I (|2007| ). However, if we 
assume a more depleted HI atmosphere with respect to 
DIV1 model, then a population of energetic HI atoms 
is required in order to fit the observations (e.g., Figures 
7 and 10). Our sensitivity study shows that depend- 
ing on the hot source position, which may be confined 
inside or outside the Roche limit, distinct spectral fea- 
tures appear at ~ ±100km/s in the wings of the transit 
absorption profile that only future high-resolution ob- 
servations with a decent S/N should be able to detect 
(e.g., Figure 7). Among the two potential sources of hot 
HI atoms that may explain the Ly-a observations dur- 
ing transit, the one confined near or outside the plane- 
tary Roche lobe is of particular interest because it may 
be directly related to the stellar wind. Indeed, if this 
hot HI population (region III) is of stellar origin, its pa- 
rameters may help constrain the stellar wind speed and 
flux. In most Hl-depleted models considered here, the 
required external source has an HI column in the range 
of [HI]hot ~ (2 — 4) x 10 13 cm~ 2 and an effective tem- 
perature in the range of Thi/Tb ~ 90 — 110 (see Figure 
10). 

In summary, our sensitivity study shows that multi- 
species observations including HI and heavy elements are 
the key tools essential for understanding the structure 
and composition of transiting exoplanets and their in- 
teraction with the nearby parent star. With all predic- 
tions expounded here, we look forward to the upcom- 
ing HST/Cosmic Origins Spectrometer (HST/COS) and 
STIS observations that should finally reveal the true bal- 
ance between thermal and non-thermal populations and, 
we hope, their spatial distribution. 
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